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Abstract 

We set up the effective field tfieories wfiicfi describe tfie S'0(5)-invariant picture of the 
high-Tc cuprates in various regimes. We use these to get quantitative conclusions con- 
cerning the size of 5'(9(5)-breaking effects. We consider two applications in detail: (z) the 
thermodynamic free energy, which describe the phase diagram and critical behaviour, and 
(n) the Lagrangian governing the interactions of the pseudo-Goldstone bosons with each 
other and with the electron quasiparticles deep within the ordered phases. We use these 
effective theories to obtain predictions for the critical behaviour near the possible bicritical 
point and the pseudo-Goldstone boson dispersion relations, as well as some preliminary 
results concerning their contribution to response functions. We systematically identify 
which predictions are independent of the microscopic details of the underlying electron 
dynamics, and which depend on more model-dependent assumptions. 



1. Introduction and Summary 



One of the remarkable features of high-Tc cuprate materials is the basic connection 
they display between antiferromagnetism (AF) and superconductivity (SC). Indeed, any 
convincing theory of these materials must explain this fundamental property. 

S.-C. Zhang [1] has recently argued for a new perspective on this AF-SC connection 
within the cuprates. He identifies an approximate SO{5) symmetry of the one-band Hub- 
bard model which is believed to describe the dynamics of the pairing electrons within the 
Cu — O planes of these materials. This 50(5) symmetry contains as subgroups the S0{3) 
symmetry of spin rotations (which is spontaneously broken in the AF phase) and the 
electromagnetic SO{2) invariance (whose breaking defines the SC phase). From Zhang's 
vantage point both ordered phases arise once SO (5) is spontaneously broken, and the 
competition between antiferromagnetism and superconductivity becomes a 'vacuum align- 
ment' problem, in which the direction of the order parameter is fixed by small eff'ects which 
explicitly break the approximate SO {5) symmetry. 

Evidence for the validity of this picture comes from the observation within the SC 
phase [2] of a collective mode, centred near momentum (7r/a,7r/a), which couples to the 
spin-flip channel in neutron-scattering experiments. Its energy gap depends on doping in 
the same way as does Tc itself, taking the value £^ ~ 41 meV at optimal doping. This state 
is understood in Zhang's framework as a pseudo-Goldstone boson (pGB) [3], whose gap is 
kept small compared to the underlying electron-electron interaction energy, J ~ 0.1 eV, 
by the S0{5) symmetry. Such a state is argued to be an approximate eigenstate of the 
Hubbard model Hamiltonian in refs. [4] . 

As is true for any physical system [5], the implications of the SO (5) picture at energies 
much lower than its intrinsic scale, J, may be efficiently encoded in terms of a low-energy 
effective theory. The effective theory which does so will depend on the symmetries and 
degrees of freedom which arise in the low-energy regime, and so can differ depending on 
the regimes of temperature or doping which are of interest. 

Of particular interest are those low-energy predictions which depend only on the 
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symmetries and degrees of freedom of the low-energy theory, and not, say, on the values of 
the eflFective low-energy coupling constants. This is because such predictions are robust, 
in the sense that they are independent of the details of how these symmetries are realized 
by the underlying electron dynamics. This robustness obviously raises the stakes of the 
comparison with experiment, since disagreements cannot be attributed to small changes 
in microscopic details. 

Our goal in the present paper is to pursue the ideas of ref. [1] by systematically 
exploring the low-energy effective theories which describe it in several regimes. A short 
summary of some of our results has been given in ref. [6]. We pause here to list some of 
our conclusions, before going into more detailed discussions in subsequent sections. 

• (1) The Phase Diagram: Our first application is to the system's phase diagram. In §2, 
we consider the system's free energy, and use it to reiterate Zhang's description of the 
qualitative features of the phase diagram in the temperature-doping (T — x) plane. We 
then extend this reasoning to argue that the size of the phases in this plane are related to 
the size of the <S'0(5)-breaking interactions, and so to the gaps for the pseudo-Goldstone 
states. 

More precisely, suppose the system Hamiltonian (or Lagrangian) density has the form 
Ti- — Hinv + e^sb, where Hinv is SO{5) invariant, and e <^ 1 quantifies the magnitude 
of the explicit <S'0(5)-breaking interactions at zero doping. We argue below that the gap 
for the pseudo-Goldstone state in the AF phase is £af = 0(e^/^J), while that in the SC 
phase is Bsc = 0{e^/^J). We learn the size of e from the observed 41 meV gap of the 
SC phase: e ~ 0.01. Both the Neel temperature at zero doping (T,v), and the SC critical 
temperature at optimal doping (Tc), are predicted to to be of order e^/^J, leading to the 
(correct) expectation that these should be of order lOOK. 

Similar statements hold for the doping required at zero temperature to destroy the 
AF order and enter the SC regime. This transition is predicted for dopings of order 
Xc — e^/^ ~ 10%. In addition, a potential 'mixed' phase (MX), which is both AF and SC 
in nature, can arise between the purely SC and AF phases, again over a range of dopings 
which is of order e^/^ ~ 10%. The superconducting phase itself is expected to extend out 
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to dopings for which the SO (5) invariance is no fonger a good approximation. 

Although the prefactors to these estimates are model dependent, and so will 

vary from material to material, the general shape of the resulting phase diagram should 
be shared by all high-Tg systems, and seems to agree reasonably well with the typical 
experimental phase diagram [7]. 

• (2) Critical Behaviour: The high-Tg systems differ from traditional superconductors in 
that they display critical behaviour within a few degrees of Tc- This critical behaviour 
is the topic of §2.3. The 5'0(5)-invariant scheme predicts the usual scaling behaviour for 
the transitions from the SC or AF regimes into the disordered phase. It makes different 
predictions for the vicinity of a potential bicritical (or tetracritical) point which is expected 
should the SC and AF phases coexist for some temperatures and dopings. Zhang has 
argued that one is attracted here to an 5'0(5)-invariant fixed point, based on the d = 2 + e 
expansion for d-dimensional nonlinear sigma-models. We point out here an alternative 
possibility where the stable fixed point weakly breaks SO (5), based on the d = A — e 
expansion, and compute the critical exponents in this picture. 

• (3) Pseudo-Goldstone Boson Spectrum: Next, §3 applies SO{b) to the underlying Hamil- 
tonian for energies, temperatures and dopings corresponding to the ordered phases, and 
obtains the properties of the pseudo-Goldstone modes. As was emphasized elsewhere [6], 
since there are more features to these gap spectra than there are parameters in the general 
effective theory, it is possible to predict low-energy relations amongst the various features 
of the dispersion relations for these states. 

Model-independent, weak-coupling conclusions do not appear to be possible for the 
disordered phase, although inferences regarding this phase may be made by making more 
specific use of additional assumptions concerning the microscopic dynamics. 

Our pGB Lagrangian differs in detail from that proposed by Zhang, since it involves 
a few more terms. Although they completely agree for qualitative purposes, we argue the 
inclusion of the most general interactions is required to permit the robust inference of the 
quantitative size of symmetry-breaking effects. 
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• (4) Response Functions: Finally, §4 uses the effective theory to make some preliminary 
points concerning the pseudo-Goldstone boson contributions to electromagnetic and spin 
response functions within the AF and SC phases. In particular, the origin within the 
effective theory, of the pseudo-Goldstone pole in the SC-phase spin response is identified. 



2. Implications for Thermodynamics 

In terms of the underlying electrons, the ordered AF and SC phases of the high-Tc 
systems are distinguished by nonzero values for the order parameters ('0p(icr^)'0_p) and 
(4'1+Q^""4'p^ J where p is an electron momentum, Q = ^) and denotes the Pauli 
matrices. Long-wavelength variations of these quantities can therefore be described by the 
fields: 

riQik) (xipp+k{ia2)i^-p and oc V'J+Q+fe^"V'p, (1) 

where k is much smaller than either p or Q. 

In Zhang's framework, in addition to the SO (3) x SO (2) symmetry of spin and elec- 
tromagnetic gauge transformations there are approximate symmetries which rotate Uq and 
Us, into one another. To represent this symmetry it is convenient to group Uq and Ug into 
a real, five-dimensional quantity 

n=f!°V (2) 
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on which the extended symmetry acts by matrix multiplication on the left by an orthogonal 
five-by-five matrix: n — > O n, where O is an arbitrary five-by-five orthogonal matrix. 
This identifies the electromagnetic S0{2) and spin S0{3) subgroups as the block-diagonal 
combinations: 

SO{3) J ■ ^'^^ 
We denote by Tq, the hermitian, antisymmetric five- by- five matrices which generate SO (5), 
with the special cases of the generators of SO{2) and SO{3) represented by: 

Q = g(^' o), and f=(° X (4) 



Here q = 2 is the electric charge of the order parameter in units of the electron charge and 
t are a basis of three-by-three generators of SO (3). 

The resulting SO (5) symmetry is simultaneously subject to two kinds of breaking. On 
one hand it is broken explicitly (but weakly) to the exact S0{3) x S0{2) subgroup by the 
electron Hamiltonian, and on the other hand it is spontaneously broken to SO{A). What 
follows in this section spells out our assumptions concerning how this symmetry manifests 
itself in thermodynamic functions. 



2.1) Thermodynamic Potentials 



Given this picture we may write the free-energy density, /, as the sum of an 50(5)- 
invariant term, /inv, and a small S'0(5)-breaking term, /gb- For a slowly varying order 
parameter we may use a derivative expansion of /, leading to the following most general 
expressions for time- independent Uq and Us'- ^ 

/5/inv = ^^inv + w'^^^ (y^anQ ■ V Va^s • VaU^ u\^^ (uq ■ V "h Us ■ Van^) 

-I- W^'' (^VcHq ■ VcHq -I- VcHs ■ VcHs j -I- U^j^^ (uq ■ VcUq + Us ■ VcTls j H , 

Pfsh = Vsh + wj"'* (VaTT-Q • Va^Q - Va^s • VaUs^ + My"'' (uq ■ Va^Q - 77,^ • Va^s^ 
-I- W^"'* (VcTIq ■ VcTIq - VcTls ■ Vc^s^ + W^"'* (hq ■ VcTIq - Us ■ Vc^s^ 
-I- tt||°'' Uq ■ VaTlQ Us ■ VaTls + ttj°^ Uq ■ V c^q Us • V c^s H • 

(5) 

Here (3 = l/kT, and the ellipses represent terms involving more derivatives than two. 
The coefRcient functions, Vinv:VshT'UJ^_i\w^^\u^^''^^ and ^y^''^'* are potentially arbitrary 
functions of the 5*0(3) x 150(2) invariants Uq ■ Uq and ■ ng. They depend as well on 
the two thermodynamic variables: temperature and doping (more about this dependence 
below). Keeping in mind the anisotropy of the cuprate systems, separate coefficients are 



This expression differs from that of refs. [1], [8] in three ways. Besides working in the isotropic hmit, 
these authors: (1) impose the constraint ns-ns+nQ-nQ=l, (2) keep only three of the lowest terms in an 
expansion of the w's and u's in powers of fields, and (3) keep only quadratic terms in the potential v. (Their 
quartic interactions arise from the chemical-potential dependence of the terms involving two derivatives of 
^ UQ.) We comment further on these differences as they arise at subsequent points in the text. 
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included for derivatives, a = x,y, parallel to the copper oxide planes, and those, c = z, 
perpendicular to this plane. For simplicity we assume rotational symmetry within the 
planes. 

In the limit of SO {5) invariance, Vsh = = w^^^ = u^^^ = ■Uy"-' = tt^^ = u^^^ = 0, 
and each of the remaining functions VinY , Wj^\ w^^\ u^^"^ and u^^^^ depend only on the 
combination n ■ n = Uq • Uq + Ug ■ ris- 

The predictive power of the approximate symmetry comes from computing observables 
perturbatively in the small symmetry-breaking interactions. Doing so requires a quanti- 
tative characterization of the size of the symmetry-breaking terms in /. We imagine, 
therefore, the microscopic electronic Hamiltonian to have the generic form: 

T-C = Hinv + e 'Hsb - A* Q, (6) 

where Hinv preserves 50(5), while Tish and /j, Q both break it. /j, here is the chemical 
potential for the electric charge, Q. We take = to correspond to half filling. 

Notice that SO (5) symmetry is explicitly broken in two diflFerent ways in eq. (6), and 
so there are two different parameters, e and /x, which govern the symmetry-breaking terms 
in the free energy, /: 

• Zero Doping: Things are simplest at half filling, for which n = and so only one 5*0(5)- 
breaking parameter exists: e. Based on our later discussion of the spectrum of the gap in 
the pseudoGoldstone boson spectrum in the SC phase, we take in what follows e ~ 0.01. 
It is this small parameter which ultimately determines the size of the different phases in 
the system's phase diagram. 

• Nonzero Doping: The second source of explicit 5'0(5) breaking in eq. (6) is the chemical 
potential, /i, which is introduced to describe the doping of high-Tg systems away from half 
filling. This breaks 5*0(5) because n couples only to electric charge, which is just one of 
the ten 5*0(5) generators. For thermodynamic applications it is more convenient to trace 
the dependence of the free energy on the number density of charge carriers, n, rather than 
fi, and so we define the doping, x,hy n = x/V, with V the volume of a unit cell. In typical 

7 



high-Tc systems at zero temperature, x ^ 10% is characteristic of the AF phase, while 
0.1 < X < 0.4 is representative of the SC phase [7]. 

To explore the implications of these two sources of symmetry breaking we consider 
first the classical phase transitions which are implied by the choice of an SO{5) invariant 
order parameter. 



2.2) Classical Phase Transitions 



For a classical phase transition the free energy density, f{x,T,nQ,ns), is assumed to 
admit a Taylor expansion in powers of Uq and Ug near Uq = Ug = 0. This assumption 
of analyticity at zero field is known to fail in the vicinity of a critical point, where the 
long-distance fluctuations can cause / to acquire a singular dependence on its arguments. 
For temperatures closer to Tc than a few degrees, the resulting critical behaviour is poorly 
described by a classical discussion. Unlike traditional superconductors, the correlation 
length in high-Tg systems is sufficiently small to permit critical behaviour to be observed, 
although typically only within a few degrees of [9] . 

For these reasons we expect a classical analysis to be adequate for the purposes of 
identifying the overall size occupied by the various phases within the phase diagram in the 
T — X plane. We investigate the implications of SO{5) for the critical behaviour in §2.3, 
below. 

We start by rederiving Zhang's treatment of classical phase transitions. Our addi- 
tion to his argument is the quantifying of the size of the symmetry-breaking effects in 
this analysis. Consider therefore expanding the potential v to quartic order in the order 
parameters, Uq and n^: ^ 



V = t'oo+^(^^20 nQ-nQ+VQ2 ns-ns^+^ V40 (ng-ng)^ -1-^22 nQ-UQUs-Us+vo^ {us-Usf 



+ ■■■, 
(7) 



This potential corrects a minor error [10] in the potential in the original version of ref. [6] , although 

2 

none of the results of this reference are altered by this change. 
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where the coefficients, Vjk, are regarded as functions of the thermodynamic variables, x 
and T, as well as the small symmetry-breaking parameter e. Particle-hole symmetry ^ 
further implies vjk must also satisfy Vjk{—x) — Vjk{x), if T and e are held fixed. 

Next divide this expression into its 5'0(5)-invariant and -breaking parts, v = v-mv+Vah, 
where: 



P 



X 



Vinv = ^^00 + 2 {J^Q ■ riQ + Us ■ Us j + - ■ Uq + Us ■ Us J H , 



(8) 



The parameters in these expressions are related as follows: vo2 = P ~ 9-, '^20 = P + 9^ 
f4o — X + k + h, fo4 — X + k — h and V22 — X — k. 



Phase 



Minimum 



Conditions 



Normal (N) 




= 


V20 > 0, V02 > 


Superconductor (SC) 


IuqI 7^ 0, \ns 


= 


V20 < 0, V02 > 


Antiferromagnet (AF) 


IuqI = 0, \ns 




V20 > 0, V02 < 


Mixed (MX) 


\ns\ 7^ 0, \nQ 




^^20 < 0, ^02 < 0, /O 


Table I 



The possible minima of the free energy density, v, each distinguished by the sign of ^02 and V20- Two 
possibihties arise if both of these parameters are negative, depending on the sign taken by the quartic 
couphng, k. fc>0, defines a minimum for which both \ns\ and \nQ \ are nonzero. If fc<0, then there is both 
AF and SC type minima, separated by an energy barrier. The corresponding phase diagram in the p—g 
plane is given for both choices for k in Fig. (1). 



Minimization of this potential gives four kinds of phases depending on whether 
and/or Uq vanish at the minimum. These four alternatives are controlled, in first approxi- 
mation, by the signs of the coefficients of the two quadratic terms, V20 and V02, as outlined 
in Table I. The case where both V20 and V02 are negative, further subdivides into two 



Alternatively, the same conclusion is also drawn in later sections by considering the free energy to 
be a function of 1^ rather than x, and using the antisymmetry of Q. See also [6]. 
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alternatives which are distinguished by the sign of the quartic couphng, k. If A; > then 
both Itt-qI and \ns\ are nonzero at the minimum, while if < then two types of minima 
coexist, which differ according to whether it is \ns\ or Itt-qI which is nonzero. When k < 
these two minima are separated by an energy barrier. 




k>0 k<0 
Figure 1 

The phase diagram for the free energy function described in the text in the p — Q plane, where p 
is the 5'0(5)-symmetric, and g is the S'0(5)-breaking, quadratic couplings in the free energy. The thin 
lines represent second-order transition lines, while the fat line is a first-order phase boundary. The two 
possibilities portrayed are distinguished by the sign of the quartic coupling, k, in the free energy. 

The corresponding phase diagram, in the p — g plane, is drawn in Fig. (1). When > 
there are four phases, and all phase boundaries define second-order transitions. The four 
phase boundaries intersect at a tetracritical point, which is defined by p = ^ = 0. When 
A; < there are only three phases, with the antiferromagnetic (AF) and superconducting 
(SC) regions meeting at a first-order transition. This phase boundary intersects the second- 
order boundaries with the disordered phase (N) at a bicritical point, again at p = = 0. 
The picture which emerges is the semiclassical phase diagram described in ref. [1]. 

More information becomes available, however, once the dependence of the transition 
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regions on the small symmetry-breaking quantity e ~ 0.01 is included. To see this, imagine 
examining v in the vicinity of half filling, where x <^\. Keeping in mind that all <S'0(5)- 
breaking interactions must vanish if both x and e do, and that Vjk must be even functions 
of x, we write in this limit: 

g = e go+g2 x^+gA • • , h = e /io+^2 x^+h4 x^+- • • , k = e ko+k2 x^+k4 x'^+- ■ • , 

(9) 

where the coefficients of this expansion are functions of T which can be unsuppressed by 
additional powers of e or x. 

Since the S'0(5)-invariant couplings need not vanish with either e or x, for them we 
instead write: 

X = Xo + X^x^ + X4x'^ + ---. (10) 

The 5'0(5)-invariant quadratic term must be handled more carefully, however. This is 
because real high-Tg systems are antiferromagnets at zero temperature and doping, and so 
we have the important additional information V2o{T = a; = 0) > and vo2(T' = a; = 0) < 0. 
Equivalently this implies: g{T = x = 0) > and -g{T = x = 0) < p{T = x = 0) < g{T = 
a; = 0), which is only consistent with g{T = x = 0) — 0{e) if p{T = a; = 0) is also 0(e). 
We take, therefore, at zero temperature: 

p = epo + P2x'^ + P4x'^ -{ . (11) 

• Transition Dopings: 

These expansions determine the domains of doping over which the various phases are 
possible, so long as x is small. For temperatures below the Neel temperature - and so for 
which —go<po<go^ define the doping above which the AF order is lost, x^^p 

is then determined by the requirement that vo2 = p — ^ = 0, for which eqs. (10) and (11) 
imply: 

xl, = -e(>^)+0(e^). (12) 
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So long as both numerator and denominator are both order unity, and p2 > g2, we see 
the prediction is: x^p = 0(e^/^) ~ 10%, which agrees well with what is observed for real 
high-Tc systems. 

Similarly, a superconducting phase arises for dopings greater than Xsc, defined as the 
value of X for which V20 = p + g = 0. Again eqs. (10) and (11) give: 

^L = -e f^)+0(a. (13) 
\P2 + g2j 

and so Xsc is also predicted to be 0(e^/^) ~ 10%, so long as numerator and denominator 
are 0(1) and P2 + 5'2 < 0. 

If we instead work to 0{x ) in V20, then a second root develop beyond which 

superconductivity is again lost. As is easy to see, this second root arises for x ~ 0(1), and 
so lies outside of the domain of the small-a; expansion, and for x potentially large enough 
to invalidate the approximate SO{b) symmetry. Notice that this implies that the optimal 
doping, a^opt — defined as the doping for which Tg is largest — is also 0(1), and so is 
unsuppressed by powers of e. 

Finally, if Xap > Xsc, and if /c > 0, then dopings satisfying Xsc < x < Xaf give both 
a superconductor and an antiferromagnet (as in the MX phase) at zero temperature. This 
phase can extend over a range of dopings which is at most as large as Xaf—Xsc = 0(e^/^). 
Otherwise, if A; < 0, the AF to SC transition is first order. 

• Transition Temperatures: Similar arguments may be applied to estimate the size of the 
critical temperatures for the various ordered phases if the T dependence of the coefficients 
in V is known. In mean-field theory the large-temperature limit of the temperature de- 
pendence of the free energy may be argued on dimensional grounds. For large T in d 
dimensions, the quadratic coefficients, p and g, are proportional to T^, with a prefactor 
that can be 0(1) for p, but which is suppressed by e or a;^ for g. 

For sufficiently large T it follows that eventually p > \g\, implying a transition to the 
disordered phase (N) . The transition temperature as a function of doping may be estimated 
by asking when the contribution is the same size as is the zero-temperature values of 
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for p and g. Since the zero-temperature limits of both p and g are suppressed by powers 
of e and/or x^, the transition temperature, T^, is predicted to be smaller than might be 
expected based on the underlying electronic scales, J ~ 0.1 eV ~ 1000 -ftT. At zero doping 
only e controls the size of symmetry breaking, so the Neel temperature, T,v, is expected to 
be 

= 0(e^/V) ~ lOOX. (14) 

Once again this is the right order of magnitude for the cuprates. Of course, these esti- 
mates can be further reduced by other effects. For example, very anisotropic systems are 
effectively two dimensional, in which case 'w^^\ u^'^'^^ — > 0, implying Tc 0. 



T T T 




Figure 2 

The phase diagram for the free energy function described in the text in the temperature-doping plane. 
As for Fig. (1), the thin lines represent second-order transition lines, while the fat line is a first-order phase 
boundary. The two possibilities for which the AF and SC phases coexist are distinguished by the sign of 
the quartic coupling, /c, in the free energy. 

We are thus led to the phase diagram of Fig. (2) in the temperature-doping plane, 
which reproduces well the generic features of the phase diagram for real high-Tg systems. 
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2.3) Critical Behaviour 



Since the symmetry of the order parameter has important imphcations for the critical 
exponents, and since some critical behaviour is experimentally accessible for the cuprates, 
we pause to explore critical phenomena more carefully here. 

As is well known, in d < 4 dimensions the Gaussian fixed point is unstable in the 
infrared, and so the critical exponents are controlled by some other infrared fixed point 
of the renormalization group. In the absence of more reliable methods, one is reduced 
to expanding about either the upper {d = 4) or lower {d = 2) critical dimensions. Our 
approach here is to expand in powers of (i — 4, and we do so for both the critical lines 
separating the disordered (N) phase from the others, as well as for the bi- or tetra-critical 
point should these lines meet. We imagine in what follows that all five of the fields in Ug 
and Uq are free to fluctuate. Our treatment of the tetracritical point therefore differs from 
Zhang's, since he draws his conclusions for this point using the nonlinear sigma model, for 
which Us ■ Us + Uq ■ Uq = 1, based on an expansion about d = 2 [11]. 

Along the critical line the S0{5) picture implies the usual 0(A'")-invariant infrared 
fixed points, where iV = 3 for the Neel transition, and N = 2 (XY universality [12]) for 
the superconducting transition ^ . Both are covered by the usual result to leading order 
in the expansion in e = 4 — d: 



where denotes the fixed-point value of v^q for the N/SC transition (or of fo4 for the 
N/AF transition). ^ u is the standard critical exponent for the correlation length, which 
diverges like ^ ~ as t = T/T^ — 1^0. There is indeed good experimental evidence [9] 
that u for the N/SC transition, is given as predicted for XY universality. 

The fixed point behaviour in the bicritical (or tetracritical) region differs qualitatively 
from that found along the two critical lines. This is because the 0( A?") -invariant fixed 



Due to the relevance of the electromagnetic coupling, the SC transition crosses over from XY univer- 
^ sality to that of a charged fluid, although undetectably close to [12]. 

Very near the critical point (and so in this section only) we ignore the system's anisotropy and rescale 
^ all fields to canonically normalize the derivative terms in the free energy. 
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point is stable only for N < 4 — 0(e), and so does not apply for AT = 5. The fixed point 
appropriate for A?" = 5 may be found numerically, and the results obtained agree where 
they can be compared with previous treatments [13], as well as with the approximate 
analytic expressions given below. 

In order to identify the fixed points it is useful to consider a model for which the 
field Us spans an A^g-dimensional space, and Uq spans a space of Nq dimensions. The 
free energy is assumed to be 0{Ns) x 0{Nq) invariant. In this case useful approximate 
expressions for the fixed points may be derived in an expansion in powers of n = A^^ — Nq 
over N = Ns + Nq. In the case of present interest Ng = S and Nq = 2, and so the 
parameter n/N = 1/5 is reasonably small. 

The stable fixed points which emerge at leading order in e = A — d then are: (i) 
the 0(A'")-invariant fixed point, for N < A; (ii) the 'decoupled' fixed point (for which 

{V22)* = 0), for {Ns + 2){Nq + 2) > 36; and (Hi) a third fixed point for intermediate values 
of Ns and Nq. It is this third fixed point which is appropriate for Ng = 3 and Nq = 2, 
and it gives: 



~.8A„, ^^^l^i^Z^X,, i^«2(7V-4)A„, (16) 



4ir2 4ir2 JV-2 4,7r^ 



2 

where Xq = e/{N + 32). For N = 5, n = 1 and d = 3 these expressions agree well with the 
values we obtained numerically to linear order in e: A*/(47r^) ~ 0.141, h^/ (47r^) ^ 0.00573, 
and A;*/(47r^) ^ 0.0340. Notice the hierarchy /i* ^ /c* <^ A* which is satisfied at this fixed 
point, and is consistent with approximate SO{5) invariance, albeit with some symmetry- 
breaking couplings which are larger than 0(1%) of their 5'0(5)-invariant counterparts. 

Since the picture of how the critical lines merge near the bicritical point differs some- 
what from that of ref. [1], we next present this in more detail. The running of the two 
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quadratic couplings, p and in the far infrared is given near the fixed point by: 



= -2p + I 2(iV + 2)A* -n/i* p+ (iV + 4)/i, - 2nA*] 



-2p + ^^^^ { [24iV(iV -2)-n\N- A) 
1 



A 



I (TV + 4)/i* - 2nA;*] p + AX^ + 2{N + 2)k^ - nK 
= -2g + {-Sn(N - 4) ^ + [aN(N - 2)^ - n^N - 4)] g]. 



(17) 



The principal directions of the flow, and the corresponding scaling exponents, near the 
fixed point are given to good approximation by: 



/ 3n(A^-4)2 \ / n[(N- 

p+ n— =^ \ q p_ q — \ 

^ \AN{N -2){^- N)) ^ ^ ^ V4iV(^- 



2)(8- iV) 



St p (18) 



and: 



1 ^ 6iV 

— ~ 2 - ^2 ^) 

^+ iV +32 



TV +32 



(19) 



For N = 5 and n = 1 these become pj^ c:i p — 0.017 and p^ g + 0.22 p, while u+ ~ 0.68 
and zv- ~ 0.58. 

The largest of these exponents defines the temperature exponent, z/ = z/_|_ ~ 0.68. The 
scaling form for the free energy density, / = t~'^'^ T (f^-), defines the crossover exponent, 
0, which is given hy cj) = vj^/v-. We have: 



N + 27V + 64 

= 

2(iV -3iV + 32) 



1.2. 



(20) 



n\{N-S)^-AS] \ „_ 22o 
AN(N-2){S-N)) P-U-^^P- 



Since > 1, the two critical lines approach the bicritical point tangent to one another, 
and to the line p_ = — or, equivalently, to the line g ■ 

It would obviously be of great interest to examine experimentally the critical behaviour 
at the bicritical point, since this would more decisively distinguish the symmetries of the 
order parameter at this point. Unfortunately, to our knowledge, there is currently no 
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experimental information concerning this critical behaviour, likely due to the disorder 
which is introduced in this region by the doping process. 

3. The Nonlinear Sigma Model: Pseudo-Goldstone Modes 

More microscopically, another regime for which 5'0(5)-invariance provides definite 
quantitative predictions is deep within the ordered phases. Here the low-energy dynamics is 
governed by the Goldstone, and pseudo-Goldstone modes, together with their interactions 
with the other low-energy degrees of freedom. Much can be said about these low-energy 
properties because low-energy Goldstone-boson dynamics is largely determined purely by 
the pattern of spontaneous symmetry breaking [5], [14], [15]. 

The present section is devoted to developing the quantitative description of this 
regime. Section §3.1 first describes the most general low-energy effective Lagrangian for 
this phase, which is then used in §3.2 to determine the dispersion relations for the pseudo- 
Goldstone bosons themselves. Their interactions are used in §4 to draw preliminary con- 
clusions about their contributions to response functions in the AF and SC phases, which 
also involves their interactions with the low-energy electron- like quasiparticles. 

3.1) The Effective Lagrangian 

Imagine, then, that T and x are chosen to lie within one of the ordered phases, and 
integrating out all degrees of freedom which involve energies larger than some scale, A. 
We wish to write down the effective Lagrangian governing the degrees of freedom which 
remain at still lower energies. Although this effective Lagrangian cannot yet be derived 
from the microscopic physics, we are guaranteed to capture its physics so long as we use 
the most general possible Lagrangian involving the low-energy modes and respecting all of 
the symmetries [5]. 

When a global symmetry, G, is spontaneously broken to a subgroup, H, the self- 
couplings and spectrum of the resulting (pseudo-) Goldstone states is described at low 
energies by the nonlinear sigma model for the quotient space G/H [5], [14], [15]. When 
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G = S0{5) and H = 50(4) this implies the lowest terms in the derivative expansion of 
the Lagrangian for this system are completely determined by two constants. Of course, 
more possibilities arise once explicit 5'0(5)-breaking interactions are introduced. 

The most general such Lagrangian involving two or fewer derivatives is again built 
from the fields Us and Uq, but with the important difference, relative to §2.1, that these 
now satisfy the constraint ns-ng + nQ-riQ = 1, since we are interested in only the Goldstone 
and pseudo- Goldstone modes. In the absence of doping, the result takes a form which is 
similar to eq. (5), C = Ci^Y + >Cgb, with: ^ 



11 
2 



J2 , , 



Csh = -V + ft A dtUQ ■ dtUQ + B dtUs ■ dtUs + C {n^ ■ dtrif^) 

- ff\ D^VariQ ■ VaTiQ + Vans • Vans + F\\{nQ ■ Vang)^ 



(21) 



+ 



where ft, f\\ and f± are constants, while V^, Ay _|_, j_, Cy _l, _l, i^H j_ and are 
potentially arbitrary functions of the S'0(3) x SO{2) invariants Uq • Hq and Us ■ ris- They 
also can depend on the temperature, T, since this can appear in C through the process of 
integrating out the high-energy modes. 

Fewer terms appear in eq. (21) than in eq. (5) because of the constraint ■ + 
■ Hq = 1 which is enforced in (21), but not in (5). Should the 'length' of the SO{5)- 
breaking order parameter also describe a propagating mode in the low-energy system, as 
might be appropriate near the critical point, then this constraint may be relaxed. Since this 
mode is not guaranteed to be in the low-energy theory far away from the phase boundaries 
— unlike the pseudo-Goldstone states — we do not include it further in this section. 

Far from the critical region the correlation length perpendicular to the planes becomes 
smaller than half the interplane spacing, and we take the system to be approximately two- 



Our notation here follows that of ref. [6], generalized to the anisotropic case. Other couplings, such 
^ as e'^'^''drnns-(dnfisxdins), are also possible in specific dimensions (in this case d=2). 
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dimensional, corresponding to the limit ^ /_l = 0. We imagine working in this limit in 
what follows, and so from here on drop the redundant subscript '||' from the coefficient 
functions -D, E and F. 

Couplings to long-wavelength electromagnetic fields is incorporated in the effective 
theory through the usual substitution daUQ — > {da — ieAaQ)nQ. (A special case of this 
coupling is the dependence on the chemical potential, /U, which enters C through the re- 
placement 9tnQ — > {dt—ienQ)nQ [17], [18].) For distances shorter than the electromagnetic 
screening length, a = eqft, this electromagnetic coupling includes the Coulomb interac- 
tions of the pseudo-Goldstone bosons. It is important to realize, however, that the strong 
short-ranged Coulomb interactions need not be included in this way, since these are inte- 
grated out to arrive at C in the first place. Although strong microscopic interactions such 
as these would complicate the derivation of the C from first principles, they play no role 
when using C at low energies. 

An important consequence of these observations now follows. When using the La- 
grangian of eq. (21), the key observation is that all of the interactions are guaranteed to be 
weak at low energies, justifying a perturbative treatment. This is because all interaction 
terms are suppressed by either a derivative or a small 5'0(5)-breaking parameter or both. 
In particular, the pseudo-Goldstone spectrum may be obtained from C in mean field theory 
by expanding in fiuctuations about minima of the potential, V. 

For later purposes a useful parameterization (which identically solves the constraint 
Uq ■ Uq + Us ■ Us = 1) is given by polar coordinates on the four-sphere: 



(although care is required to properly handle those points where these coordinates are 
singular). In terms of these variables, and including a chemical potential, the Lagrangian 

A more detailed modelling of the dimensional crossover within the context of Bose-Einstein conden- 
^ sat ion is given in ref. [16]. 





(22) 
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becomes: 



f2 



(l + 2A sin^ d + 2B cos^ 9 + 2C sin^ d cos^ ^) {dtdf 
+ (1 + 2A) cos^ 6'(af(/) + eqiif + (1 + 2S) sin^ ^(^(ata)^ + sin^ a{dt(3f^ 

- (i + 2D sin^ ^ + 2E cos^ ^ + 2F sin^ 9 cos^ ^) ( Va^)^ 

+ (1 + 2D) cos^ ^(Va0)^ + (1 + 2E) sin^ ^ ((VaQ;)^ + sin^ a{VaPf)] +■■■, 

(23) 

where all coefficient functions are to be regarded as functions of cos^ 9. 

To this point we have not yet used much information concerning the nature or size of 
the explicit symmetry breaking. This we now do by making an assumption as to how the 
symmetry-breaking terms transform under SO (5). Since there are two types of symmetry 
breaking, a choice must be made for each. 

• Doping: Doping has been incorporated into the effective theory through the chemical 
potential, /i. We return to the connection between /j, and x in §3.3 below. 

• Intrinsic Breaking: Symmetry breaking also occurs at half filling, with strength e. The 
resulting symmetry-breaking pattern is S0{5) — > S0{3) x 50(2). We assume this to be 
done with the simplest possible 'order parameter', M, which transforms in the adjoint 
representation of SO{5). That is, choose M — ediag (3, 3, —2, —2, —2). 

The Lagrangian is then the most general function of the fields n = (^^), A*(5 ^^nd M, 
subject to the following SO (5) transformation property 

£{On, OfiQO^, OMO^) = £{n, /iQ, M), (24) 

where O is an SO{5) transformation. 

The utility of identifying Q and M may be seen when C is expanded in powers of 
the small quantities e and /i. Since M and Q always appear premultiplied by these small 
numbers, this expansion restricts the kinds of symmetry breaking which can arise order by 
order, which in turn constrains the possible ^-dependence of the coefficient functions in C. 
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For example, a term in the scalar potential involving 2n powers of n must have the 
following form: 

^in)= "' XI Ck:,h,...,k^,l 

(25) 

Only even powers of Q enter here due to its antisymmetry, and the term ki = li = is 
excluded from the sums due to the constraint n^n = 1. Expanding £ to low order in the 
5*0 (5) -breaking parameters e and /x necessarily also implies keeping only the lowest powers 
of Uq ■ Uq = cos^ 9 in V. 

Similar conclusions may be obtained for the other coefficient functions in the La- 
grangian of eq. (21). Working to 0(e^, e//^,//"^) in V, and to 0{t,ii^) in the two-derivative 
terms then gives: 

V = Vq + V2 cos^ + )-Va cos^ 9, (26) 

and 

A = + ^2 cos^ 6', S = So + ^2 cos^ 6', C = Co, 

(27) 

D = Dq^D2Qjos^9, E = Eq + D2Cos'^ 9, F = Fq, 
for the coefficient functions in eq. (21). Notice that the terms proportional to cos^ 9 in A 
and B are identical, as are the corresponding terms in D and E. Expanding in powers of e 

and fj,, the constants in eqs. (26) and (27) start off linear in e and //^: Ai = Aj^e+A^^ iJ,'^-\ 

etc.. The only exceptions to this statement are: Bq, Eq oc e (no fi^ term), Co, Fq oc /x^ (no e 
term), and V4 = Vi^e'^+Vi^e/i'^+V^^ii^. Furthermore, since the fi^ Uq -Uq term in V arises 
from substituting dt ^ dt — iefiQ in the kinetic term for Uq, we have: 1^2°"'^ = ~|/t^e^g^ 
to leading order. Higher powers of fi originate from terms in C which involve more than 
two derivatives. 

In this way we arrive at an effective lagrangian very similar to that of ref. [1] (in the 
isotropic limit). The main difference here is the power counting of the symmetry break- 
ing terms. Ref. [1] keeps terms quadratic in Uq and has three susceptibility parameters 
controlling the time-derivative terms, and three stiffness constants governing the spatial 
derivatives. Although a quadratic scalar potential captures the leading order in e, it does 
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n^(eM)'=i(//Q)^'in ■■■ n^(eM)'=" (//(5)^'"n 



not appear that a three-parameter derivative term corresponds to any fixed order in e or 



3.2) Pseudo-Goldstone Dispersion Relations 

We now turn to the calculation of the pseudo-Goldstone boson dispersion relations. 
The scalar potential of eq. (23) has three types of extrema: 

(1) ^0 = or tt; 

(2) ^0 = - or — ; 

(3) ^0 where c = cos^o satisfies V'{c^) = 0. 

This leads to the four classical ordered phases found in §2.2: (i) SC phase: extremum 
(1) is a minimum, and (2) is a maximum; (ii) AF phase: (2) is a minimum, and (1) is a 
maximum; (in) MX phase: both (1) and (2) are maxima, and (3) is a minimum; or (iv) 
metastable phase: both (1) and (2) are minima, and (3) is a maximum. This analysis 
becomes identical to that of §2.2 if V is assumed to be quartic in cos 6, as was done in 
ref. [6]. 

• Superconducting Phase: An expansion about the superconducting mimimum, = 0, 
gives the dispersion relations in this phase for the four bosons. The result is a spin triplet 
of pseudo-Goldstone modes for which 



m = 



c^ + 



(28) 



with the phase speed, c^gb{'^^)i S^P; ^Igb{^^) = ^Ic given to lowest order in 
S'0(5)-breaking parameters by: 



c%ASC)Jj,[l + 2{E{l)-B{l)) 



fi 



l + 2(Eo-Bo] +2(D2-A 



2 _-2V'il) 



(29) 



^sc 



-2(^2 + Vd 
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In both of these results the first equation uses the general effective theory, eq. (23), while 
the second equality incorporates the additional information of eqs. (26) and (27). 

The remaining field, (f), describes a bona fide gaplcss Goldstone mode. Its dispersion 
relation, E{k) is a more complicated function of c^/c^ and eq/j, whose form [18], is not 
required here. Its phase velocity, = c^g{SC), is given by 



clsiSC)=^ l + 2{Dil)-Ail) 
-11 



(30) 



l + 2(^Do -^o) + 2(p2-A2 



Recall in these expressions that y(cos^ 9) includes any /x-dependent contributions com- 
ing from the kinetic terms, or their higher-derivative counterparts, and a prime denotes 
differentiation with respect to cos^ 9. 

• Antiferromagnetic Phase: Expanding about the AF minimum gives the usual two magnons 
satisfying dispersion relation (28), with: 



cls{AF)^-l 1 + 2[E{0)-B{0) 
Jt 

f! 

Sl,{AF) = 0. 



1 + 2 



(31) 



The remaining two states form a pair of electrically-charged pseudo-Goldstone bosons 
satisfying: 

1 

c k + t ± eqfi, 



with: 



E±{k) = 



claB{AF) = 



(32) 



/i 



1 + 2(d(0)-^(0)) 

1 + 2(Do-^o) 

2V'{0) 

2V2 

fr 
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(33) 



These expressions imply the simple formulae of ref. [6]: = m — k/j, , and = 
-m^ + K/x^ - ^^^4^ where = 2Vi^e/f^ + 0{e'^), k = -2V§^ / + 0(e) = e^q^ + 0(e) 
ande^2]402//2^O(,). 

Within the AF phase the pseudo-Goldstone boson gap is seen to fall linearly with 
p?: fa ^af{^)[i^'af ~ A*^]) where /Iaf represents the doping for which one leaves the 
AF regime. Similarly varies quadratically with By eliminating parameters of the 
Lagrangian in favour of properties of the gap as a function of n we find the relations of 
ref. [6]: 



.2 f,,^ 4.(0) 



J2 



2 2 
l^AF-V^ 



Fopt (34) 

£L(Q) _o 4c(opt) 

2 — 2 ' 

^AF A'opt 

where //opt here denotes the chemical potential corresponding to the maximum gap, Esc- 
We expect this to occur at optimal doping, /^opt = A*(2^opt)- 

Similarly, the phase velocities for all modes in both SC and AF phases are equal to 
one another, and to ft/fp in the strict <S'0(5)-invariant limit. (The parameters ft and 
/|| are related to the compressibility and magnetic penetration depth in the next section.) 
The 0(e) corrections to this limit also satisfy some model- independent relations, which 
follow by eliminating parameters from the above expressions: 

cliSC) - cliAF) = cliSC) - cliAF) = 0(e). (35) 



3.3) The Connection Between x and fj, 

The previous expressions giving the dependence of physical quantities in terms of 
the chemical potential, /U, would be more useful if expressed in terms of the physically- 
measured quantity, the doping, x. This relation is determined implicitly in the present 
section. 
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The dependence iJ,{x) is found by adjusting /j, to ensure that the net electric charge 
equals x charge carriers per unit cell. This must be done differently in the AF and SC 
phases. 

In both phases the total charge is carried by both the charged pseudo-Goldstone 
states, and the ordinary electron-like quasiparticles responsible for conduction. The exis- 
tence of these electron-like quasiparticles at low energies in the SC phase is demanded by 
the evidence in favour of a dxi-yi-wave gap in the cuprate superconductors [19]. These 
experiments argue for the existence of ungapped states due to the nodes of the dx2-y2-waye 
gap function, which cannot be provided by the four pseudo-Goldstone modes. In the AF 
phase these degrees of freedom correspond to ordinary unpaired electrons. 

• AF Phase: In the AF phase the charge is carried by a mixture of electrons and charged 
pseudo-Goldstone states. The condition of equilibrium between these two types of charge 
carriers implies their electric chemical potentials must be equal, and so the defining con- 
dition for ij,{x) becomes: 



ex f°° ( I 1 r 1-1 

(Pem)=^ = J du ^J\fp{u) np.{ii) -np{-fi) +Nb{uj) nB{n) -nB{-l^) |, 



(36) 



where J^p and N'b respectively denote the density of states for fermions and bosons. For 
weakly interacting bosons in d spatial dimensions, whose dispersion relation is E'^ = fc^c^ + 
the density of states is given explicitly by Mb — {0,d/{2n)'^) k'^~'^E/c^, where Vtd is 
the solid angle swept out by a vector in d spatial dimensions (so VL2 — 2tt and l^s = Att). 



The Bose-Einstein distributions are given by np{ii) = 



and UBifJ') 



I ■'of small fj, it is the fermions which dominantly contribute to p, and the dependence 
x{fj,) which results is (for d = 2 space dimensions) linear: x oc /i/J. This linear dependence 
changes once /j, becomes of order the scalar gap at zero doping, £^af(0) = 0{e^/'^J), since 
at this point the scalar charge density varies strongly with signalling the transition to 
the condensed (SC) phase. 

§2.2 showed the doping for which one exits (at low temperature) from the AF phase to 



25 



be, in order of magnitude, = 0{e^/^). We see that this corresponds to a chemical 

potential whose size is /x = /Iaf = 0(e^/^J). (An identical conclusion regarding the size of 
Haf may be drawn from the condition for Bose-Einstein condensation: ~ £af{0)-) 

We see that since the chemical potential depends linearly on doping within the AF 
regime, eqs. (34) implies a linear dependence of £'ip{x) on x'^: 

SlM=£UO)i^lF-x^)- (37) 

Corrections to this linearity in x'^ arise as the SC phase is approached since /j, no longer 
varies linearly with x near the point where the bosons condense. 

• SC Phase: Deep within the SC phase we assume the charge density to be dominated 
by the Uq condensate. Writing the potential as V{n) |(m^ — i''^)n? + ■^g'^ri'^ (where 
V = eq^ft, m? = 2V2iii = 0) ^ 2V^^e, and g'^ = 12F4(// = 0) ^ 12Vi^e'^), we see (pem) is 
given by 



. , _ eqx dV 



n=no ^ 



(u^ - m^) , (38) 



where no = 6(^'^ — m^)/g^ minimizes V. This leads to the asymptotic expressions v fa 
(firVeV/t)'/^ (X a;V3 if (p^^) /(Gf^m^) > 1, and z/ - m gx/{12Vftm'') oc a; if 
g^ (Pem) /{Qftm^) « 1. 

In the case of interest we have (for d = 2 space dimensions) ft = 0(1), rn'^ = 0{eJ'^), 
and g'^ = 0(eV), and so (pem) /(Qftm^) ~ xe^/'^/{6VJ'^). Taking J ~ 0.1 eV and 
V ~ (10l)2 we find 6VJ^ ~ 10-^ leading to g'^ (pem) /{Qftm^) ^ lO^xe^^ ^ iq4^ g-j^^,^ 
this is much greater than unity when x > Xsc = 0(e^/^), it follows that /i oc x^^^ within 
the SC phase. 

We may now determine the size of fx at optimal doping. We found in §2.2 that optimal 
doping occurs for x = Xopt = 0(1), and so Xopt/xAF = 0{e~^). But since jj, oc x^^^ within 
the SC phase, we see popt ~ PAFixopt/xAF)^^^ — 0{e^/^). 

Together with the previous results, eqs. (34), these expressions imply that Sg^{x) is 
quadratic in the variable fj/^ oc a;^/^. Similarly, the gap in the SC phase at optimal doping 
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is related to the AF gap at zero doping by: 



£sc{opt) 



(39) 



AF 



This last equation, together with the interpretation of the 41 meV state as the pseudo- 
Goldstone boson of the SC phase, and the underlying electronic scale J ~ 0.1 meV, give 
the order of magnitude of the symmetry-breaking parameter: e ~ 1%. 



In order to measure the properties of these pseudo-Goldstone particles, it is necessary 
to understand how they contribute to the spin and electromagnetic response functions of 
the materials. This is the topic of the present section. Because these particles are weakly 
coupled, their response may be computed perturbatively. 

As we shall see, it shall become important for these purposes also to understand how 
the pseudo-Goldstone states couple to the other degrees of freedom in the low-energy sys- 
tem. For this reason we also write down the electron/pseudo-Goldstone particle couplings 
in this section. 

The starting point for calculating the response functions is to identify the dependence 
on the (pseudo-) Goldstone bosons of the spin and electromagnetic currents. These are 
very easily obtained, to lowest order in the derivative expansion, by constructing the 
corresponding Noether currents using the Lagrangian of eq. (21): 



Because the pseudo-Goldstone bosons are weakly coupled, correlations of these cur- 
rents may be be directly evaluated for free bosons, Uq and Ug, plus perturbative corrections. 



4. Response Functions 




(40) 
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This perturbative evaluation conveniently organizes the contributions to the response ac- 
cording to which states are responsible. The following sections give some examples of such 
calculations. 



4.1) High-Energy Contributions 

Even though the effective lagrangian only contains as degrees of freedom the states 
which actually appear in the low-energy spectrum, it nonetheless carries the information 
as to how states at higher energies contribute to response functions. The contributions of 
higher-energy states are incorporated as they are integrated out to produce the effective 
lagrangian itself, and so their contribution to correlation functions may be read directly 
from the lagrangian. 

More precisely, imagine coupling to a long- wavelength electromagnetic field, or to 
a fictitious field, s^, coupling to spin, by making the substitution: daUQ — > (da—ieQAo,) Uq 
or daUs — > {da — iTaS'^) Us, with Ta and Q the matrix generators of SO{?i) and SO{2), as 
defined in eq. (4). The high-energy contribution to the electromagnetic and spin response 
may then be obtained by differentiating the effective Lagrangian twice with respect to A^^ 
or s^. The resulting correlation is proportional to 5'^(x — x') 5{t — t'), as would be expected 
for a mode which fiuctuates on a scale much shorter than that over which the response is 
computed. 

Making this substitution for the electromagnetic response leads to the magnetic pen- 
etration depth along the planes. A, and the electric screening length (or compressibility), 
a. These are given by: 



where fig is the magnetic permeability of the material. The result in the ground state is 
obtained by using the ground-state configurations: Uq ■ Hq = 1 (SC) or Uq ■ Uq — (AF). 

Similarly, the spin response obtained in this way gives the high-energy contribution 
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to the paramagnetic susceptibility: 



(xab) =filf^{l + 2B)nsTanns, (42) 

\ / h.e. 

where fig is the magnetic moment of the pseudo-Goldstone bosons. These depend on 
temperature, only through the weak dependence of coefficients ft, f\\, A, B and D. 

4.2) Goldstone Poles 

The next simplest contributions to response functions to compute are the poles which 
occur in the correlation of the currents for spontaneously broken (approximate) sym- 
metries, due to the contribution of the corresponding (pseudo-) Goldstone states. This 
includes the superconducting contribution to the electromagnetic response in the super- 
conductor, and the magnon contribution to the spin response in the antiferromagnet, 
in addition to the pseudo-Goldstone boson couplings to the additional SO{5) generators 
which rotate the spin and charge degrees of freedom into one another. These are given in 
perturbation theory by the Feynman graph of Fig. (3), which describes the direct creation 
and destruction of the Goldstone boson from the ground state by the current of interest. 

X ^ 

Figure 3 

The Feynman diagram which produces the Goldstone pole contribution to the current-current corre- 
lation function. The dashed lines represent pGB propagation. 

The result for the SC electromagnetic and the AF spin response, obtained by summing 
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the pole and the high-energy contributions, is: 



^^(^)([Pem(x,^),Pem(0)]) 



f d?k dw 



—i(jjt+ik-x 




GB 



J (27r)2' 



2 



^^W<[Pspin(^,^),Pspin(0)]> 



GB 



f cPkdw 



~iuit+ik-x 



(1 + 2B)f^Pclf jnsTgriXrjnns) 
-uj'^ + k'^clf -i^afi^ 

(43) 



where Cgc = Cgb{SC) and Caf = Cgb(^F). and 7^/ represents small damping con- 
tributions to the pGB dispersion relations. For the spin-spin correlation, rj denotes the 
normalized eigenvector of the gap matrix corresponding to the GB direction whose pole is 
being computed. Notice eqs. (43) reproduce eqs. (41) and (42) in the limit a;,/c ^ 0. 

Of course, eqs. (43) are not specific to S'0(5)-invariant physics, since they also hold 
for any antiferromagnet or superconductor. This is because they simply contain the im- 
plications of the Goldstone bosons associated with spontaneously broken 5'0(3) or 5*0(2) 
invariance. As a result, although they provide a good description of the response in these 
phases [20], [21], [7] this is not a real test of /SO (5) symmetry. ® It is the contributions 
of the psettdo-Goldstone states in both phases which provide more interesting informa- 
tion. After pausing to address a puzzle concerning the pseudo-Goldstone pole in the spin 
response of the SC phase, we close by describing some of the features of the pGB response. 

4.3) Pseudo-Goldstone Poles: A Puzzle 

It was the contribution to neutron scattering of the spin-triplet pGB state in the SC 
phase which originally motivated the SO{5) picture. Since neutrons couple to electron 
spins, the 150(5) interpretation of the neutron-scattering experiments requires the pGB to 
contribute a resonance to the spin response functions. 

This immediately leads to a puzzle. The spin-triplet pGB contributes a pole similar to 
eqs. (43) to three 5*0(5) currents (Zhang's n operators) which are spontaneously broken in 
the SC phase. But Fig. (3) does not produce a pole in the SC spin response function in the 

For the same reasons, neither are such responses detailed tests of the explicit models in which they 

Q 

are usually derived. 
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SC phase, because the symmetry of spin rotations is not broken in this phase. The puzzle 
is how such a pseudo-Goldstone pole can arise as a resonance in the SC spin correlation 
function. This section sketches how this puzzle is resolved within the effective field theory 
framework. 

The difficulty with producing a pGB pole in the spin correlation function lies in the 
observation that pspin of eqs. (40) involves only even powers of the boson field, n^. The 
same is true of the interactions in the effective lagrangian, eq. (21), and so it is difficult 
to generate a graph of the form of Fig. (4), which would generate a resonant contribution 
to the spin correlation function. The 'blobs' of Fig. (4) represent any graphs which can 
produce the triplet state starting from quasiparticles created by the spin density, Pspin- 

In our opinion the resolution of this puzzle comes from the couplings of the pseudo- 
Goldstone bosons to the electronlike quasiparticles in the SC phase. As stated earlier, the 
existence of these quasiparticles can be inferred from the evidence for cia;2_j^2-wave pairing, 
since some of these experiments indicate the existence of gapless excitations in the SC 
phase. These excitations must be in addition to the four Goldstone and pGB states of the 
SC phase. 



Figure 4 

The Feynman diagram which produces the pole contribution to the spin correlation function in the 
SC phase due to the spin-triplet pseudo-Goldstone state. The 'blobs' represent fermion loops, while the 
dashed line represents pseudo-Goldstone boson propagation. 



To see how these quasiparticles can help with the puzzle, suppose them to have the 
quantum numbers of electrons, and to have couplings that are weak, so that their propa- 
gation is approximately described by the Lagrangian density Cq = f d'^p cjj[idt — Ep Cp. 

Here Cp = ( j destroys a quasiparticle which propagates with dispersion u> = Ep, as 
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is described by Cq. Standard arguments [22], [23] can now be used to identify those in- 
teractions which are the most important in the long- wavelength limit. Goldstone boson 
couplings are all irrelevant (in the renormalization group sense) in this limit, but the least 
irrelevant of their couplings to the electronic quasiparticles involve the emission and ab- 
sorption of a single Goldstone particle. The resulting electron Lagrangian density which 
describes this is: 

9 sip, k) cJ^fc^Cp • nQ{k) + gQ{p, k) cl_^_f,{ia2)c*_p nQ{k)] + h.c. (44) 

This interaction is least irrelevant for special regions of momentum of the quasiparticle 
pairs [22] , such as when the net momentum of the pair is close to zero. 

Notice that an expectation value for Uq introduces a quasiparticle gap, proportional 
to Qq (uq), so a (ia;2_j^2-wave symmetry of the gap restricts how gQ{p,k) can depend on 
momenta lying on the Fermi surface. Approximate SO (5) invariance relates the coupling 
5'q(P) k) to gsip, k), and so implies a similar dx2_y2-waive symmetry for gsip, k). 

The principal observation at this juncture is that these quasiparticles contribute 
quadratically to the spin density, Peispin oc J cfip c^pffcp, and so the coupling of eq. (44) 
introduces a pGB pole into the spin correlation function through the Feynman graph of 
Fig. (4), with the 'blobs' representing quasiparticle loops. Even though this graph does 
not arise at leading order in perturbation theory, the singular shape of the pole permits 
it to dominate the lower orders for energies and momenta which are related by the pGB 
dispersion relation. 

4.4) Pseudo-Goldstone Boson Response 

Away from special features, such as the poles just discussed, the dominant contribution 
made by pseudo-Goldstone bosons to response functions arises through the Feynman graph 
of Fig. (5). 
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Ant= / d'^pd'^k 



Figure 5 



The Feynman diagram which produces the pseudo-Goldstone boson contribution to the electromag- 
netic response in the AF phase. 



At zero frequency and momentum transfer this graph contributes a temperature- and 
doping-dependent contribution to the response functions, which are precisely those due to 
the noninter acting gas of bosons having the dispersion relation of the pseudo-Goldstone 
states. Since this dispersion relation is 'relativistic', the results are those of a gas of 
relativistic bosons [24] . 

In the SC phase one finds in this way the thermal paramagnetic susceptibility due to 
the spin-triplet pGBs: 



XpGb{SC) = 



dx 



X 



d-3 



E 



X 



2 + (d-2)E2] ns(iU = 0). 



(45) 



Recall is the solid angle swept out by a vector in d spatial dimensions (O2 = 27r, 



1^3 = 47r), and the boson dispersion relation is E 

-1 



sections, nB{lj) 



x"^ + S^, for X = kc. As in previous 
^f3(E-eqn) _ ^ denotes the Bose-Einstein distribution function. The 
thermal electric screening length due to the charged pGBs of the AF phase is given by a 
very similar expression: 



(AF) 



pGB 



{2ncy 



dx — — — 
E 



x'^ + {d- 2)^2 Us ifi) + ris i-fi) 



(46) 



Both Xpgb(5'C) and (l/o)pQB i^F) are therefore seen to be exponentially activated, 
oc e~'^^, for kT <^ £, and to vary as T^-^ for kT :$> £. Their contribution to the specific 
heat per unit volume, c^^, of the corresponding phases is also exponentially small for 
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kT <^ £, and varies as T"^ for kT ^ £. Unfortunately, the exponential suppression makes 
this T dependence diflBcult to detect at low temperatures, while the large-T power-law 
behaviour only applies, as derived, for T much greater than the pGB gap, and yet small 
enough that the sample remains in the ordered phase. 

The dynamic response function of a relativistic bose gas is also known for nonzero fre- 
quencies and momenta [25] . This carries considerably more information about the pseudo- 
Goldstone boson response, although the differentiation of the pGB contributions from other 
degrees of freedom is easiest at low frequencies. We close by presenting some preliminary 
remarks concerning this response, and defer a more detailed applications to experiments 
to a later publication. 

Because of its electric charge, the pGB state of the AF phase cannot contribute to a 
diagram of the form of Fig. (4), and so give a pole at low temperatures in the electromag- 
netic response. Fig. (5) nevertheless does produce some strong dependence on frequency, 
due to the singularity it implies near the threshhold for producing pairs of pGBs. Below 
this threshhold the electromagnetic response from Fig. (5) is purely real, while it is complex 
above. As a result, for temperatures T > S^f the pGBs contribute zero to the conductivity 
for frequencies below threshold, u> <^ 2Saf, but the conductivity then grows steeply beyond 
this threshold. The implications of the resulting expressions for electromagnetic scattering 
from cuprates in the AF phase will be reported elsewhere. 

4.5) The Disordered Phase 

Many of the results obtained above for pseudo-Goldstone bosons deep in the ordered 
phases might also be expected to apply in the disordered phase. If so, the wealth of 
experiments available there would permit many more detailed tests of SO (5) symme- 
try. Applications to the normal phase are also theoretically appealing, since a number of 
striking features might be expected within the SO{5) picture, including two of the more 
striking implications pointed out in ref. [1]: the explanation for the pseudogap, and of 
the connection to a successful scaling analysis of the temperature dependence of NMR 
relaxation times [26]. Furthermore, the absence of spontaneous breaking of SO{5) also 
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implies that the boson contribution to correlation functions of the SO (5) currents, j^{x) - 
with a = 1, . . . , 10 labelling the SO{5) generators - are very simply related in the S0{5)- 
invariant limit: (j^ j^) oc Sat- 
in this section we make some caveats concerning the use of the effective lagrangians of 
previous sections in the disordered phase. Our main point is to emphasize that conclusions 
drawn from the effective lagrangian involving tiq, Hs (and possibly electronic quasiparti- 
cles) are not protected in the disordered phase by the general low-energy constraints of 
Goldstone's theorem, and so are necessarily more dependent on assumptions made about 
the details of the underlying electronic interactions. Although this makes these predictions 
no longer simply consequences of the symmetry-breaking pattern, they can nevertheless 
be worthwhile as sources of more detailed information about this underlying microscopic 
physics. 

We next describe some of the ways in which model dependence can enter predictions 
made for the disordered phase using the boson lagrangians described in this paper. 

• 1. Degrees of Freedom: First, the system's real degrees of freedom need not be as 
assumed, since no general principles require the existence of low-energy bosonic states de- 
scribed by Uq and Us if no spontaneous symmetry breaking occurs. An exception might be 
in the immediate vicinity of the transition line into one of the ordered phases, since conti- 
nuity at this line would require the gap for the pseudo-Goldstone states of the disordered 
phase to still be small in the ordered phase. The deeper one moves into the disordered 
phase, the less one would generically expect the boson gaps to remain small compared to 
the intrinsic scale, J. 

Furthermore, since the four pseudo-Goldstone states should fill out a linear repre- 
sentation of the unbroken SO{5) of the disordered phase, fluctuations in the modulus of 
the five-dimensional vector, n, should also appear in the low-energy spectrum. This ar- 
gues that the effective lagrangian of interest is of the form considered here for the free 
energy, (eq. (5) supplemented by time-derivative terms) rather than as was used for the 
pseudo-Goldstone bosons of the ordered phases (eq. (21)). These differ through the relax- 
ation of the constraint Us ■ ris + Uq ■ Uq = 1, due to the fluctuations in the magnitude of 
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5'0(5)-breaking order parameter, n. 

• 2. Weak Coupling: Even if the boson degrees of freedom exist at low energy, they need 
not be weakly coupled since (unlike Goldstone bosons) they are not required to decouple 
at low energies. It is noteworthy, however, that for weakly-coupled systems the electron- 
boson interactions given in eq. (44) are the only couplings between the electrons and bosons 
which can be marginal or relevant (in the RG sense of refs. [22]). By contrast, there are 
a number of self-couplings among the bosons which can be relevant or marginal in the 
infrared. Of course, the existence of strong couplings among the low-energy degrees of 
freedom doesn't necessarily invalidate the use of the effective lagrangian, it could just 
complicate the extraction of its predictions. 

• 3. Bose-Einstein Condensation: Even if the previous assumptions should apply to a par- 
ticular system, it is still true that a weakly-coupled version of electrons and bosons cannot 
provide a good description for the cuprates in the disordered phase for dopings larger than 
optimal. This is because if the bosons Uq are supposed to appear at low energies, and if 
their couplings are weak, then the relation, iJ,{x), between chemical potential and doping 
should be reasonably well described, as in §3.3, by a gas of free bosons and electrons. But 
this description always implies Bose-Einstein condensation for sufficiently large dopings, 
since for large enough x the bosons always 'win' and, by condensing, dominate the ex- 
pression for the electric charge density. This cannot describe the observed decrease of the 
critical superconducting temperature with increasing doping, above optimal doping. 

Further work is necessary to better explore these implications for the normal phase, 
and to more clearly identify which predictions of the effective lagrangian for this phase are 
model-specific, and which are more robust consequences of the symmetry-breaking pattern. 
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